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Abstract — Over the last years, novel low-complexity approach- 
es to the equalization of MIMO channels have gained much 
attention. Thereby, methods based on lattice basis reduction 
are of special interest, as they achieve the optimum diversity 
order. In this paper, a tutorial overview on LRA equalization 
optimized according to the MMSE criterion is given. It is proven 
that applying the zero-forcing BLAST algorithm to a suitably 
augmented channel matrix — the inverse of the square root of the 
correlation matrix of the data symbols times the noise variance 
forms its lower part — results in the optimum solution. This fact 
is already widely used but lacks a formal proof. It turns out that 
it is more important to take the correlations of the data correctly 
into account than what type of lattice reduction actually is used. 
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I. Introduction 

The joint reception of signals transmitted in parallel — either 
considering multi-antenna systems or multi-user scenarios — 
will become even more important over the next years. When 
designing transmission systems for such multiple-input/mul- 
tiple-output (MIMO) channels, the interference amoung the 
individual signals has to be dealt with by means of equaliza- 
tion. 

During the last decade, numerous techniques known from 
intersymbol-interference channels — e.g., linear equalization, 
decision-feedback equalization (DFE, also known as succes- 
sive interference cancellation (SIC) and also used in the 
Bell Laboratories space-time (BLAST) approach), maximum- 
likelihood detection, cf. J5] Table E.l] — have been transferred 
to the MIMO setting. However, novel approaches based on 
lattice basis reduction, e.g., ETIl . ifTTl . are of special interest. 
Using these lattice-reduction-aided (LRA) techniques, low- 
complexity equalization achieving the optimum diversity be- 
havior lfT31 is enabled. 

In this paper, a tutorial overview on LRA equalization 
optimized according to the minimum mean-squared error 
(MMSE) criterion is given. It is shown that it is more important 
to take the inherently introduced correlations of the data 
symbols correctly into account, than which lattice reduction 
approach actually is used. The main result of the paper 
is to establish a connection of the V-BLAST algorithm to 
the successive MMSE estimation of correlated data starting 
from basic principles. It is proven that applying (zero-forcing) 
BLAST to a suitably augmented channel matrix — having the 
inverse of the square root of the correlation matrix of the data 
symbols times the noise variance as its lower part — indeed 
results in the optimum solution. To the best knowledge of the 
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authors, a formal proof for this fact yet has not being presented 
in literature. However, it has been used widely without taking 
care of its validity. 

The paper is organized as follows: in Sec. [XT] the channel 
model is introduced and conventional equalization techniques 
are briefly reviewed. Lattice-reduction-aided equalization is 
addressed in Sec. [HI] and its MMSE DFE version is analyzed 
in detail in Sec. [IV] Concluding remarks follow in Sec. [Viand 
in the Appendix fundamentals on estimation are compiled. 

II. Channel Model and Equalization 

We consider uncoded multiple-antenna transmission over fiat- 
fading channels where joint equalization at the receiver side 
is possible. The input/output relation is given by the usual 
equation^ 

y = Ha + n . (1) 

This model has either to be understood in complex base-band 
notation or as its real-valued model with doubled dimension- 
ality ]9). As all subsequent discussions can either be applied 
to the complex or the real model, we do not distinguish both 
approaches in the sequel. In each case, the channel matrix 
is expected to be of dimension Nr x Nt- The differences 
between both views are discussed if required. For successive 
schemes, due to the larger degree of freedom, usually the real- 
valued model has some advantages J4). 

Each component of a is independently drawn from a 
zero-mean one-dimensional M-ary ASK constellation A = 
{±l/2,±3/2,...,±(M - l)/2} or an M 2 -ary QAM con- 
stellation, with an M-ary ASK per quadrature component. 
The correlation matrix of the data vector hence reads <fr aa = 
E{aa H } = cr^J, with variance a 2 a = E{|a M | 2 }. The noise is 
assumed to be spatially white with variance cr 2 per component, 
i.e., $„„ = E{nn H } = er 2 I. 

A. Linear Equalization 

The interference between the parallel data streams can be 
eliminated by means of equalization, i.e., via r = H^y a 
decision vector is generated. Having r, individual threshold 
decision can be performed. 

Using Unear equalization (LE), optimized according to the 
zero-forcing (ZF) criterion, the receive matrix reads 



11 



(LE.ZF) 
R 



(h h h) 



(2) 



'Notation: A T : transpose of matrix A; A H : Hermitian (i.e., conjugate) 
transpose; A~ H : inverse of the Hermitian transpose of a square matrix A; I: 
identity matrix; matrices are denoted by uppercase letters, vectors by lower 
case letters. E{ }: expectation. 
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i.e., the receive matrix is given by the Moore-Penrose left 
pseudo inverse of H. Already in [9| it has been observed that 
the minimum mean-squared error solution is obtained by using 
the augmented matrix (£ = f % is the inverse signal-to-noise 
ratio) 



H 



H 



y = 



(3) 



(ATr+JVt) 



(N K +N T )xN T 

in the ZF solution and feeding y into the resulting receive ma- 
trix rather than y. Subsequently, all quantities corresponding 
to the augmented channel model are marked by a horizontal 
bar. 

B. Decision-Feedback Equalization 

Some gains over linear equalization can be achieved by 
using sorted decision-feedback equalization, also known as 
BLAST or SIC. The required matrices for ZF DFE are 
obtained by performing a sorted QR-type decomposition such 
that 

HP = QL , (4) 

where P is a permutation matrix (a single one in each row and 
column), Q is unitary and L is lower triangular. From these 
quantities, the feedforward matrix F and the lower triangular, 
unit main diagonal feedback matrix B are calculated as^| F = 
diag(.L) -1 Q H and B = diag(Z) -1 .L, respectively. 

Again, the MMSE solution is obtained by plugging the 
augmented channel matrix into |@}, cf. JTSJ. 

III. Lattice-Reduction-Aided Equalization 

Unfortunately, using linear equalization or DFE, only a diver- 
sity order of Nr — Nt + 1 (for the complex-valued model) 
is possible. Lattice-reduction-aided equalization schemes, e.g., 
lETl . ifTTl . have proven to require only low complexity, never- 
theless being able to achieve the full diversity order of the 
MIMO channel lfT31 . The idea is to choose a "more suited" 
representation of the lattice spanned by the columns of the 
channel matrix H; equalization is done with respect to the 
new basis, which is desired to be close to orthogonal. At the 
very end, the change of basis is reversed. 

A. Lattice-Reduction-Aided Linear Equalization 

For performing LRA equalization, in the first step lattice 
basis reduction, e.g., by using the LLL algorithm lfl2l (or 
some complex-valued version thereof, e.g., [6)), is performed 
to obtain 

H = CZ , (5) 

where Z € 2 JVtxJVt i s a n integer unimodular matrix, i.e., 
has only integer coefficient^ and |det(Z)| = 1. The reduced 
channel matrix C is usually required to have columns close 
to orthogonal and of small norms (depending on the definition 

2 diag(-) denotes a diagonal matrix with elements taken from the main 
diagonal of the indicated matrix. 

3 In case of complex signals, the set of integers Z has to be replaced by the 
set Z + jZ of Gaussian integers. 



of "reduced"). Using (0, the receive signal is given by y = 
CZa + n. 

In the second step, only C is treated and the signal 
z = Za, which is taken from a translate of the integer lattice 
(ZZ Nt = Z, Nt ) and hence can be obtained by individual 
threshold decision per component, is to be estimated. This 
transformed data vector has zero mean, fi z = E{z} = 
E{Za} = ZE{a} = ZO = 0, but is correlated with 
covariance matrix 



* 22 = E{zz H } = E{Zaa"Z H } = a 2 a ZZ H 



(6) 



Third, the change of basis is reversed via Z~ . 

1 ) LRA ZF Linear Equalization: Applying LRA ZF linear 
equalization the correlations are ignored and the receive matrix 
is simply the left pseudo inverse of the reduced channel matrix 



H 



(LRA,LE,ZF) 
R 



(7) 



2) LRA MMSE Linear Equalization: As in the conventional 
case, the MMSE solution may be obtained by applying all 
operations to the augmented channel model, cf. Il20l . Hence, in 
the first step H is fed into the lattice basis reduction, resulting 
in (note: Z usually differs from the ZF case) 



H = CZ 



(8) 



Using the definition of H, the reduced augmented matrix can 
be written as 



C = 



H ' 




HZ 1 ' 


def 


c 




z- x = 






A 



(9) 



with the obvious definitions of C and A. The receive matrix 
(with respect to y) is then given by 
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(LRA, LE, MMSE) 
R 



H^_ 1(5 H 



(C c 



(10) 



or with respect to y, when deleting the last Nt columns 



H 



LRA, LE, MMSE) 
R 



= (C H C + (Z- H Z- 1 ) C H (11) 



(12) 



Z [H U H + QI) H H 



This receive matrix takes the correlations of the data per- 
fectly into account. To see this, note that from the basic 
literature on estimation, e.g., lfT4l Theorem 2.6. 1], the optimum 
MMSE linear estimator is given by 



Vr 1 



(13) 



which, since white channel noise was assumed and <& 22 = 
ZZ H , exactly gives the receive matrix (fTTl i. The covariance 
matrix of the resulting minimum mean-squared error e is given 
by 



( c H * _1 c 

' c H c + cz~ H z~ 1 



(14) 
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y 

0=3 



Fig. 1. Lattice-reduction-aided DFE. 



fi. Lattice-Reduction-Aided ZF DFE 

In order to enhance performance, linear equalization can be 
replaced by DFE, resulting in lattice-reduction-aided DFE, cf. 
FigHJ 

As in the classical case, for performing DFE, the (sorted) 
QR-type factorization of the respective channel matrix is 
required. For LRA ZF DFE, the factorization has the form 



CP = QL 



(15) 



Feedforward and feedback matrices F and B are calculated 
as explained above. In the feedback loop, the components of 
z are detected in an optimized order described by the per- 
mutation matrix P. After reestablishing the original ordering, 
an estimate of the original data vector a is generated via the 
inverse of the integer unimodular matrix Z. 

IV. Lattice-Reduction- Aided MMSE DFE 

The optimization of the LRA DFE according to the MMSE 
criterion is not as straightforward as in the ZF case. This is 
due to the correlation of the data symbols Zk to be estimated 
in an optimum succession within the DFE loop. Up to now, 
in the literature this fact has not been treated in detail; 
usually simply the ZF solution with respect to the augmented 
matrix has been used, e.g., Il20l . lfT3l . We first review the 
straightforward application of the BLAST algorithm (8) to 
the augmented channel model and then compare these results 
to those obtained from the theory of optimum estimation of 
correlated Gaussian random variables. 

A. Lattice Reduction 

As in the LRA MMSE linear case, we stick to the aug- 
mented channel model H and consider the lattice reduction 
according to (0 and (|9). Assume for simplicity of notation, 
that the columns of C are sorted according to the optimum 
decision order, i.e., we replace C implicitly by CP, thereby 
anticipating the permutation matrix P to be determined during 
the calculation of the required matrices. Thereby, the optimiza- 
tion criterion is — as proposed in the V-BLAST system — the 
noise enhancement encountered in the feedforward processing. 
For the MMSE solution this criterion is identical to looking at 
the minimum main diagonal element of the error covariance 
matrix. 

B. V-BLAST Algorithm 

We first simply perform the (MMSE) V-BLAST algorithm 
with respect to the augmented channel matrix C. Assuming 
that I (I = 0, . . . , Nt — 1) symbols are already known, the 



BLAST approach is to simply delete the I first columns (due 
to the assumed sorting) of C and proceed with the residuafl 
augmented channel matrix Cm. 

1 ) Feedforward Matrix: Having deleted the first I columns, 
the potential feedforward matrix (with respect to the aug- 
mented channel model) for estimating the remaining Nt — I 
symbols reads 



=,(') 



<0 



nil) 

JNt-1 
C [l] C V\ ' 



C [l] C [t] 



A V\ A W 



-1 r 



r H 4 H 

C [2] A [l] 



(16) 



In each step the row f^', corresponding to the symbol Zk 
which can be detected most reliably, is appended to the entire 
feedforward matrix F. The feedforward matrix for the non- 
augmented, original channel is obtained from F by deleting 
the last Nt columns. 

2) Optimum Sorting: In the BLAST algorithm, usually the 
norms of the row of the feedforward matrix are considered 
as sorting criterion (8). These are proportional to the noise 
enhancement and hence determine the error rate. Using ( TT6l >. 
these row norms are given by the diagonal elements of 



f «\f (1) ) h 



c [i]C[i] 

c [i] c m 



c 



H ( 



c [i] c [i] 



(17) 



If C has already been sorted optimally, the upper left 
diagonal element will be the smallest. Otherwise, the first 
row of F and that with the smallest norm are exchanged; 
this exchange is also recorded in the permutation matrix P. 
After Nt iterations the entire feedforward matrix F and the 
optimum processing order, represented by the permutation 
matrix P are known. 

3) Feedback Matrix: Knowing F and P, the feedback 
matrix B can be calculated. It is well-known [7| that the 
approaches of a) canceling before applying the feedforward 
matrix (as usually proposed in the BLAST context) and b) 
canceling at the output of the feedforward matrix (as is pre- 
ferred in the DFE context) are equivalent. Here, we consider 
the latter strategy, cf. also Fig. [TJ 

Since it is optimum to cancel all known interference, the 
feedback matrix calculates to 

bi 



B = 



'Nt- 



= FCP = F 



(18) 



As in each step FmCm = I holds, i.e., the remaining symbols 
are equalized and the already canceled are ignored, it is easy to 
see that B is a lower triangular matrix with unit main diagonal. 

4 Given a matrix M, let A<frn denote the matrix obtained from M by 
deleting the first I columns, and A^rn denote the matrix composed of the 
first I columns of M, i.e., M = [JVfrri Please distinguish that from 

the notation ■">, which indicates a quantity present in step I (counting from 
zeros to Nt — 1) of the algorithm. 
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Moreover, by construction, the rows of F are orthogonal; via 
a diagonal gain matrix G we can write F = GQ , where Q 



The conditioned covariance matrix <J>~ 



can be written as 



is an (TVr + Nt) X iVx matrix with orthonormal columns. In the one hand 
summary, using the lower triangular matrix L = G~ 1 B < fT~8T > 
can be written in the form 

CP = QL , (19) 



deP' Zl 

follows. Since from © = A = [Ax A 2 ], we have on 



1 



—t;A A 

(7 



i.e., applying the BLAST algorithm a sorted QR-type (QL) 
factorization of the reduced augmented channel matrix C is 
inherently performed. 

In more detail, the (I + l)th row of the feedback filter is 
given by 

bi+i 



fx CP , 



(20) 



which, using (1161 1. is the first row of the matrix 



M {1) = 



Writing CP 
can write 



C[t\C {l] 



A\ An 



A "i A m 



with Co 



C m A {1] 



C[j] and A 2 



P 

An 



(21) 



C 2 C*2 



A 2 A 2 



C H 2 A H 2 



C\ c 2 
A x A 2 



(22) 



C 2 C 2 



A 2 A 2 



C" 2 C X + A^A, 



C. Optimum Estimation of Correlated Data 

We now turn to the situation of deriving the required 
matrices directly from the theory of minimum mean-squared 
estimation and the properties of correlated random vectors 
when parts of the variables are already known. Looking at the 
optimum linear estimator (l37l i. summarized in the Appendix, 
feedforward and feedback matrices can immediately be given 
by identifying the respective quantities suitably. 

However, from (T37T > it can be deduced that the optimal 
processing depends on the mean and covariance matrix of the 
vector of not yet detected symbols. These quantities, however, 
depend on the previous decisions when performing DFE. In 
turn, optimum filtering and the optimum processing order 
potentially may depend on the actual decisions made so far 
within the DFE. In the following we show, that this is actually 
not the case. All required matrices can be calculated in advance 
and the influence of previous decisions is taken into account 
via the feedback matrix in an optimum way. 

1) Feedforward Matrix: Again assume that the first I 
symbols (contained in the vector Z\) have already been 
detected. We can partition the vector z, mean vector and 
correlation matrix of this vector in the form 

*11 *12 



Z = 



Z 2 



Ml 
fJ-2 



(23) 



Under white noise, the feedforward matrix — filtering the re- 
ceive vector y (non-augmented model) for obtaining estimates 
of the remaining Nt — I symbols — is given by (cf. (|37}) 



if 



J Nt- 
C "] C il] 



C [i]C [lVi 



u n^ 22 \ Zl 



22|zi 



-,H 1 



(24) 



A»Ax 
A H A 

A 2 A X 



A U A 

H 

A 2 A 2 



(25) 



On the other hand, with the partitioning (1231 and using IflOl 
Page 472, Eq. (7.7.5)], we can write (elements marked by * 
are irrelevant) 



*21*ll 1 *l 2 ) 



(26) 



A comparison of d25T l and d26b reveals that for all I, we have 
< (*22 - *21*n*l ^ 



A^A 2 =A^A {1] . (27) 



■ 22 ^21 ^11 ^12 

Hence, the optimum feedforward matrix calculates to 

-1 



Ml) 



c [i]C m 



A "] A ll] 



-['] 



(28) 



2) Optimum Sorting: According to the general theory of 
estimation (Eqs. (l38l and (HTli). given Zi and applying the 
optimum linear estimator (feedforward processing), the corre- 
lation matrix of the error with respect to the remaining, not 
yet known symbols Zk is given as 



°i (c" 



[l]C [t] + ° 2 n {3> 22 - *21*ll 1 *l 2 ) 1 



A m A u] 



(29) 



The next symbol to be detected is the one, for which the 
corresponding main diagonal element of <I> ee is minimum. 
Assuming the channel matrix has been accordingly rearranged, 
the upper left main diagonal element is the smallest and only 
the first row of the feedforward matrix is used to produce a 
decision symbol. Otherwise, the respective rows are exchanged 
which is kept track of in the permutation matrix P. 

3) Feedback Matrix: From d37l i and using d40b . d4~Tb . the 
feedback filter follows immediately, too. The influence of the 
already detected symbols has additionally to be canceled from 
the receive vector y. This is done by remodulating the vector 
Zx of decisions via Crn, containing the first / columns of 
C. Moreover, the mean of z 2 given Z\ has to be taken into 
account (starting from fi z = 0). The task of the feedback filter 
is hence twofold: to cancel the known interference and at the 
same time to predict the not yet decided symbols from the 
known ones. 

With the goal to have the cancellation point at the output 
of the feedforward matrix, the feedback filter, when already I 
symbols are known, calculates to 



-1 

22|zi 



(30) 



C\qC {l] 



n^22|z : 



c [i] c [i] 
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and with the above abbreviations (partitioning of C), after 
straightforward manipulations, we arrive at 

= ( (c 2 i c 2 + a^aX 1 c%c 2 - iWfcu 1 



+ (C 2 H C 2 + A 2 H A 2 J C^C 1 

= (C 2 H C 2 +A 2 H A 2 ) _1 

■{C%C X -A%A 2 * 21 *^ 

From d25l l. the correlation matrix is given as 



(31) 



*P ZZ 



*11 *12 
*21 *22 



A-^A^ A^A 2 
A 2 Ai A 2 A 2 



-1 

def 


o 


V~ 


-1 




v H 


u 





o - vu^v" 



Again using iflOl Eq. (7.7.5)], we have 
$ii 

*2i = [v"o- L v - u 

and together with A 2 A 2 = U, we arrive at 



v"o~ l 



(32) 

(33) 
(34) 



A" 2 A 2 * 21 <t>^ 



u ( v H o W - u 



v H o 1 o VU~ l V H 



= u ( v H o~ 1 v - u 
'u 



(35) 



In summary, the feedback matrix, when I symbols are already 
known, is given by 

M {1) = (C^C 2 + A^A^ 1 (C%C 1 + A 2 A^\ . (36) 

Assuming that the symbols Zk are in the optimum ordering, 
as for the feedforward matrix, since only a single next symbol 
(the currently best) is decided, only the first row of the matrix 
is actually used. Note that the respective row of the 
feedback matrix B is obtained from that row by appending a 
single one and then At — I — 1 trailing zeros. 

D. Comparison and Discussion 

From the above derivations it is immediate that both per- 
spectives on LRA MMSE DFE lead to the same result. A 
comparison of ( fTSI l — here deleting the last Ar columns to 
return from the augmented to the original channel model — 
and (l28l l reveals that for both cases the feedforward matrices 
are identical. 

The sorting is based on dTvT > and 1291 . respectively. As the 
feedforward processing is identical, this also holds for the error 
variances or the norms of the filter vectors, proportional to 
these variances and hence the same decision orders result. 

Finally, the feedback filters are also identical; this is re- 
vealed by comparing d22l and ([36}. 

Hence, the straightforward application of the V-BLAST al- 
gorithm for sorted QR decomposition to the extended channel 



matrix indeed results in the optimum solution to LRA MMSE 
DFE. The "trick" behind this lies in the lower part of the 
augmented matrix. Whereas for classical DFE the (scaled) 
identity matrix is present, in case of LRA the inverse of 
the square root of the correlation matrix of the vector z 
to be estimated is present (cf. (0 and d25ll). As shown, 
deleting columns and calculating the feedforward matrix on 
this reduced channel matrix has the same effect as updating 
the correlation matrix of the residual symbols. 

The above derivation also reveals that in case of MMSE 
DFE for correlated symbols the feedback matrix fulfills two 
tasks: the cancellation of the interference of already detected 
symbols and some kind of prediction of the still unknown 
information symbols from the symbols up to now known. In 
case of white data symbols, only cancellation is required. 

Numerical simulations reveal that it is more important to 
take the correlations of the data symbols correctly into account 
than using a specific type of lattice reduction. Conducting 
lattice reduction on the original channel matrix, and using 
the resulting matrices C and Z to create an augmented 
matrix on which the QR decomposition is done, performs only 
marginally worse than starting rightaway with the augmented 
matrix. However, using the LLL on the original channel 
requires less complexity and is independent of the current 
SNR. 

As the above derivation is valid for any channel model 
and any correlation of the data, we can conclude that when 
performing MMSE DFE for correlated symbols, optimum 
feedforward and feedback matrices and the optimum sorting 
can be calculated via the V-BLAST algorithm. Thereby, the 
algorithm simply has to work on an augmented channel matrix, 
which has the inverse of the square root of the correlation 
matrix of the data symbols times the noise variance as its lower 
part. In other words, all required matrices are obtained by 
performing a sorted QR-type decomposition of this augmented 
channel matrix. However, with regard to computational com- 
plexity this procedure is far from optimum as the algorithm has 
to work on a matrix of approximately doubled number of rows. 
Fortunately, the efficient "fast V-BLAST algorithm" proposed 
in 0J Table II] can simply be modified to take correlated 
data (correlation introduced via a matrix Z) into account. 
Here, only the computation of R and Q according to 0] 
Eqs. (26) and (28)] has to be modified. Using the initializations 
(notation from Q) R = aZ H Z and Q = (l/a)(Z H Z) -1 , 
this algorithm efficiently delivers the same results as the ZF 
BLAST algorithm applied to the augmented channel matrix. 

LRA equalization for MIMO channels can be viewed as the 
counterpart to partial -response signaling (PRS) ATI . J2] for in- 
tersymbol interference channels, see J5]- In both cases an inte- 
ger polynomial/matrix is split from the actual channel transfer 
function/matrix and only the residual system is considered. 
Equalization is done towards the target polynomial/unimodular 
matrix. The non-whiteness of the data sequence to be detected 
has to be taken into account for MMSE equalization of PRS 
(e.g., J2] Appendix A]); the same is true in LRA schemes. 
However, in contrast to PRS, which is usually employed to 
achieve some desired transmitter side characteristics (spectral 
zeros at DC or Nyquist frequency), the use of LRA enables 
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full diversity of the MIMO transmission system and hence is 
the key to significantly improve error performance of uncoded 
transmission. 

V. Summary and Conclusions 

Lattice-reduction-aided equalization optimized according to 
the MMSE criterion of MIMO channels has been studied. For 
the first time it has been proven that applying the zero-forcing 
BLAST algorithm to a suitably augmented channel matrix — 
having the inverse of the square root of the correlation matrix 
of the data symbols times the noise variance as its lower part — 
indeed results in the optimum solution. It is more important 
to take the correlations of the data correctly into account than 
what specific type of lattice reduction actually is used. 

Finally it should be noted that taking the uplink/downlink 
duality |[T6l into account, instead of employing receiver-side 
equalization, MMSE LRA precoding |[T9l can be performed. 
The given results can immediately be transferred to this 
transmitter-side technique, which is of great importance in the 
multi-user downlink. 

Appendix 

Some Fundamentals of Estimation Theory 

In this appendix, for convenience, two important properties 
on minimum mean-squared error estimation of correlated 
and non-zero mean random variables are reviewed from the 
literature. 

First, we consider a vector x with (possibly) non-zero 
mean /j, x and covariance matrix & xx . This vector is observed 
through the matrix H and disturbed by (zero-mean) Gaussian 
noise n with covariance matrix <!>„„. Hence, the observation 
y = Hx + n is present. The optimum linear estimator for 
this setting is given by, e.g., [14, Page 68] 

x = (H H ^H + ^"V^-^y - H/iJ + n x 



( 



,y 



(37) 



The covariance matrix of the resulting estimation error can be 
written as 



H n <i> H 



(38) 



Second, assume a multivariate Gaussian distribution (ran- 
dom vector w) of dimension Q with mean fi w and covariance 
matrix & ww . Let the random vector, the mean vector, and the 
covariance matrix be partitioned according to 



w = 



Wl 

w 2 



Ml 

M 2 



*11 *12 
*21 *22 



(39) 



where the dimensions of the upper and left parts are q, e.g., 
dim(tt)i) = dim(/x 1 ) = q, dim($n) = q x q, etc. 

Having already knowledge on the first q components of 
the random vector w — i.e., the vector u;i — the mean and the 



covariance matrix for the residual Q — q variables (vector w 2 ), 
conditioned on the knowledge w-y, calculate to 

»2\ Wl = M 2 + *21*U (™1 - Ml) ( 4 °) 
&22\ Wl = *22 - *21*il *12 ■ (41) 

Note that the new covariance matrix is the Schur complement 
of <I>ii in <I>; it does not depend on the actual value of Wy. 
Note additionally that both quantities can be obtained in one 
step or successively in q steps, each time assuming additional 
knowledge of a single symbol. 
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